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朱 法 集 ， 鲍 媛 媛 
(江苏 大 学 理学 院 ， 镇 江 212013) 


搞 要 : 本 文 研究 了 顾客 批量 到 达 且 服务 台 忙 时 与 闲 时 故障 率 不 同 的 多 服务 台 可 修 排队 系统 ， 推 导出 
了 系统 的 稳 态 平衡 方程 以 及 稳 态 概率 值 的 求解 思路 。 由 于 N > 2 时 手工 计算 的 复杂 性 ， 使 
用 Mathematica 软件 编程 实现 了 稳 态 概率 值 的 求 取 ， 并 最 终 得 到 了 系统 有 效 服务 台数 的 稳 态 分 
布 、 稳 态 队长 的 母 函 数 以 及 平均 队长 等 重要 的 系统 指标 。 
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1 引言 


近 几 年 ， 不 同形 态 的 可 修 排队 系统 已 经 得 到 了 深入 地 研究 ， 获 得 了 很 多 的 成 果 。 关 于 多 服 
务 台 可 修 系统 的 文献 主要 有 : 文献 [1 研究 了 M/MIN 并 联 可 修 排队 系统 ;文献 2] 研究 了 有 
一 个 修理 工 ， 且 忙 时 与 闲 时 故障 率 相 同 的 M/MIN 可 修 排队 系统 ， 给 出 了 多 服务 台 可 修 系统 
稳 态 分 布 存 在 的 条 件 ; 文献 [3] 和 文献 | 各 研究 了 有 一 个 修理 工 ， 且 人 忙 时 与 闲 时 故障 率 不 相同 
的 M/MI/N 可 修 排队 系统 ， 并 给 出 了 这 类 多 服务 台 可 修 系统 稳 态 分 布 存在 条 件 的 一 种 新 形式 。 
可 以 看 出 ， 关 于 故障 率 可 变 的 多 服务 台 可 修 排队 系统 的 研究 正在 起 步 ， 而 对 批量 到 达 的 多 服务 
台 可 修 排队 系统 的 研究 尚未 见 诸 于 文献 文献 中 已 有 的 研究 仅 给 出 了 多 服务 台 可 修 排队 系统 稳 
态 概率 值 简要 的 求 取 过 程 以 及 N = 1 时 稳 态 概率 值 的 解析 表达 式 ， 由 于 N > 2 时 系统 稳 态 概 
率 值 手工 求 取 的 复杂 性 ， 没 有 给 出 其 相应 的 解析 表达 式 ， 没 有 为 相关 的 一 些 实际 问题 提供 一 个 
解决 问题 的 工具 ， 鉴 于 此 本 文 在 讨论 了 批量 到 达 且 故障 率 可 变 的 多 服务 台 可 修 排队 系统 的 基础 
上 ， 将 通过 使 用 Mathematica 软件 编程 实现 N > 2 时 稳 态 概率 值 的 求解 过 程 ， 并 给 出 了 一 个 具 
体 实 例 ， 为 现实 生活 中 这 一 多 服务 台 可 修 系 统 的 研究 提供 一 个 简便 的 工具 。 


2 模型 描述 

本 文 的 基本 假设 如 下 : 

1) 系统 中 及 个 服务 台 ， 一 个 修理 工 。 顾 客 成 批 到 达 ， 每 批 到 达 的 人 数 X 是 一 个 随机 变 
量 ， 批 到 达 流 是 泊 松 流 ， 参 数 为 和 ，X 的 分 布 为 


Oo 
P{X=i}j=7, i=1,2,.…, 》 关 = 
i=1 
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Oo 
E(X) < co， D(X) < co， R(z)=) mix 2 
+ 


同一 批 到 达 的 顾客 按 随机 顺序 接受 服务 。 服 务 时 间 服 从 负 指 数 分 布 ， 均 值 为 1/y， 到 达 过 程 与 
服务 过 程 是 相互 独立 的 。 假 设 一 服务 台 在 空闲 时 的 故障 率 为 C1,， 在 忙 时 的 故障 率 为 Cc。， 且 故障 
状态 的 服务 台 被 修复 的 概率 为 6。 其 中 入 , 1, C2, 5 和 jy 都 是 非 负 常数 。 

2) 若 服务 台 故 障 发 生 在 修理 工 空闲 时 ， 服 务 台 立 即 被 修理 ， 否 则 等 待 ， 修 理 完 毕 ， 服 务 
台 立 即 恢复 服务 ， 且 服务 台 修复 如 新 。 服 务 规则 为 RCFS， 服 务 中 断 时 ， 正 在 被 服务 的 顾客 可 
转移 到 其 他 空闲 的 有 效 服务 台 接受 服务 ， 若 没有 空闲 的 有 效 服务 台 ， 则 在 原 服务 台 处 等 待 。 

以 X() 表 示 时 刻 t 系统 中 的 有 效 服 务 人 台数，Y(t) 表示 时 刻 t 系统 中 的 顾客 数 ， 则 随机 过 
程 {X 的 ( 的 ;二 > 0 描述 了 系统 时 刻 上 的 瞬时 状态 ， 这 一 随机 过 程 是 二 维 Markov 过程。 
以 Pij(t) 表示 系统 在 时 刻 t 的 瞬时 状态 概率 ， 即 


的 = P{X(t) =i, Y(t) = i=0,b.,N; j=1,2,.…, 
用 已; 表示 系统 处 于 状态 (i,j) 的 稳 态 概 率 ， 即 


i POS Ys TON Jel er 


3 ”系统 的 稳 态 平衡 方程 


通过 分 析 这 一 模型 的 各 个 状态 的 具体 情况 ， 则 可 以 画 出 这 一 系统 的 状态 转移 图 。 根 据 状态 
转移 图 ， 我 们 可 以 写 出 系统 所 满足 的 稳 态 平衡 方程 如 下 : 
1) 当 i=0 时 
(A+0Po=GPo，7=0， 
a 


(+OP,j= GP;+ > NPprik, 了 > 0. 
k=0 


2) 当 0<i<N 时 


(i +AtPo=6P_ 10+HP1+(it+laPr, j=0, 
[jC +(i—)t+A+t+6 十 了 已, 


了 7 一 1 


=6P-17+(GO+DHURH+[G +G+T1- 力 人 ] 忆 + 二 >» APixrTjy—k, 7 < 
二 
| 
lic + A+6+ip]P,; = 6P_ 13 +ipPsrit+ [liz + G1]Pir + 》， MPipTj-k, j=%, 
k=0 


j-1 
[iG2 + A+6+ip] Ps = 6P 1 +ipPsn t+ (二 TCR 十 > APi,pTj-k, J > 
k=0 
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3) 当 ;i= N 时 
(NG + NPio = P10+ HPn, j=0, 
j-1 
[((N =— DG +jc2+Atijin]P,; = 6P_ + (7+ 1uPyrit+t 》>， NPiprj-k, 0< 了 < 也 
k=0 
3 
[Ne +A+NJPy = 6Pi1j + NpPisrit+ YO APker-e 了 > N. 
k=0 


4 ” 稳 态 概率 值 的 编程 实现 
设 N 
Gi(z) = 2Pj， G(z)= ,Gi(z), 0<i<N, lz<gl, 
;二 0 i=0 


则 G(z) 即 为 系统 队长 的 母 函 数 ，Gi(1) = 并 > Py 即 为 系统 中 有 i 个 服务 台 处 于 有 效 状态 
的 稳 态 概率 。 运 用 与 文献 罗 类 似 的 方法 ， 最 终 可 将 系统 有 效 服务 台数 和 队长 母 函 数 的 求解 
问题 转化 为 Bio, i = 1,2,… ,N， 这 N 个 概率 值得 求解 问题 ， 而 这 N 个 概率 值 可 通过 N 个 
关于 Pio, i = 1,2,… ,N 的 相互 独立 的 方程 求 得 四 ， 从 而 最 终 得 到 系统 的 性 能 指标 。 由 
于 N > 2 时 ， 手 工 计算 这 些 稳 态 概率 值 较为 复杂 ， 已 有 的 文献 仅 给 出 了 N = 1 时 系统 稳 态 概率 
值 的 表达 式 ， 对 N > 2 时 仅 给 出 了 解 的 存在 性 证 明 以 及 求解 的 过 程 ， 并 没有 实现 过 。 本 文 在 已 
有 文献 的 基础 上 ， 通 过 使 用 Mathmatica 软件 编程 实现 了 N > 2 时 本 文 所 讨论 的 这 一 类 可 修 系 
统 稳 态 概率 值 的 求解 过 程 ， 该 程序 是 使 用 Mathmatica 软件 的 基本 语句 独立 编写 而 成 ， 在 以 往 
的 文献 中 未 出 现 过 ， 主 要 程序 见 下 : 


n—l 
R= >, riz'; mt= Table[Switch [i—j, 1,—6z, 0], i, n+ i,j, n+ 1]; 


i1 
For [i=1,i<n,i=it+l, mtl,]= (i Dit+(i~ le+At+6— AR)z— (i— Dal; 
mtltr+1,n+1]= (n+nc t+ A—AR)z— ny; 

For[i=1,i<n,i=i+l, mtlli,it+1] = —ic2z]; 

b= Table[0,1,n+1]; blll,1]] = (G — C2)P,o0z; 


i—l 
For E =2,i<ni=itl bd]= 5 (mC tp G0) -PP yi misi™ 


m=1 


i a1 
十 ba m(C1 — Ca)Pii_mz "+t! 一 >», | 
?9 一 1 


m=1 


blllynt+1])= 5 (mC tp 0) -PPon-mz "td, mpPrn_-mz ™ 
m=1 


m=1 
b= Transpose[b]; For[i=1,i<nt+l,i=it+l,p=mt; pl[Al,d] = bl[AU, a]; y= 8]; 
g= NlSolvelhh == 0,z]]; f=z/.9; WU= Lengthlf; j=1; 
For6E = 1,i<U,i=it+l, HI0< fl) <1, = 加; 
Print[kj]; 2 = kj; d; = Detlyi] == 0; Print[ld;]; j = j+1]] 


h= Table [(iG + 和 +6)Pio == jpPs1 + 6Pi-10 + (i+ oP,o,i,2,n— 1]; 
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kk = (nO+ NPno == pjpPn1+éPn-10; UU= Append[h, kkl; 


For[i=1,i<n—2,i=it+l, b= (it2+At+ih+ (no—id)Ps 


i—l1 
== (i+ 1)pPriti+é6Pn_lit+ > PrpMri-k; U= Appendlll, bi 
k=0 


Forli=3,i<n—1,i=i+l], Forlj=1,i<i-2,j;=ji+1), 


cj = (j2 +AtIp+ nod)P,s == (7 +1)pPisyr + 6Pi-1; 
了 -1 


+》 Poem-z+(iG + (it+1— I )Prs; l= Appendlll, ei); 
k=0 


hh = Detlmt]; hk= Factorlhh]; pk= hk/(z—1); a1= Detlyi]; bi = Factor[ai]; 
pk2=b1/(z—1); z=1; kg = 6Go0— 62G1== (0 一 人)Poi 
itl 


hh = Table [ic， (itl)GGitt == 》 ml(G — Ga)Piriini mi bn— 1|; 


m=1 


万 一 >》， Gi==1; lls= Append[hh,kg]; lis = Appendllls,tt]}; si1={); 
i 二 0 


Forli=0,i<n,i=it+!, s1= Appendl[si,Gi]]; sst = Solvellls,s1]; sl1= s1/.sst; 
Forli=0,i<n,i=i+!1, sg = sil[1,i+ 1]]; hkl = pk* sgo == pk2; ll= Appendlll,hk1]; 
For[i=1,i<n—1,i=it+], l= Appendlll, di]]; result = Solvelii!] 


为 了 检验 程序 的 准确 性 ， 将 N = 1 代入 程序 ， 此 时 系统 中 仅 有 一 个 未 知 变量 ,1， 通 过 程 


序 计算 的 结果 为 
_ 64— 6A— Xa 


me 
得 到 的 这 一 结果 与 手工 计算 的 结果 一 致 ， 从 而 验证 了 程序 的 准确 性 。 由 于 服务 台数 较 多 以 及 该 
系统 包含 的 自 变量 众多 ， 因 而 得 到 的 包含 这 些 变量 的 稳 态 概率 值 的 解析 表达 式 非 常 巨大 ， 在 此 
就 不 再 一 一 显示 这 些 包含 变量 的 稳 态 概 率 值 的 解析 表达 式 ， 仅 就 一 具体 实例 进行 计算 并 显示 最 
终结 果 。 令 N=4 入 =2,6=4, k= 8, 1 = 1, C= 2， 将 其 代入 到 程序 中 ， 得 出 结果 如 下 


》 


Pi,o = 0.0945623, Pz.o = 0.182642， Pi = 0.0488143, 
Ps,0 = 0.230792, Ps,1 = 0.0571713， Ps,2 = 0.00828753, 


Pio = 0.222297， Pa1 = 0.0513266, Piz = 0.00603873， Ps = 0.000380259. 


5 ”系统 的 性 能 指标 


根据 以 上 求 出 的 概率 值 ， 可 求 得 Gi(1), i = 0,1,… ,N， 进 而 可 求 得 系统 的 稳 态 可 用 
度 4 = 学 人 Gi(1): 系统 的 稳 态 故 障 频 度 为 


N co 
B= ,2 ,PlG -DG + jCal; 


i=0 j=0 
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系统 的 稳 态 更 新 频 度 为 2 
CO = 》、 6Gi(1). 
i=0 


同 理 可 求 得 Gi(z), i = 0,1,.… ,N， 进 一 步 可 得 到 系统 队长 的 概率 母 函 数 


N 
G(z) = > Gi(2), 


i=0 
再 由 母 函 数 的 性 质 ， 系 统统 计 平 衡 下 的 平均 队长 为 
dG 

BlY]= Ee z=1 


另外 ， 其 它 的 性 能 指标 有 : 所 有 服务 台 有 效 时 系统 中 顾客 数 的 期 望 值 Ly = 3 jPN,;; 
有 i 个 服务 台 有 效 时 系统 中 顾客 数 的 期 望 值 Zi = 21jPi,;; 系统 中 顾客 人 数 的 期 望 值 


N oo 
Ls = ,2 iPn;; 


i=0 j=1 
系统 中 忙碌 的 服务 台数 的 期 望 值 
N o0 
E(B)= 》 >》 Pi 
i=0 j=1 
空 亲 的 服务 台数 的 期 望 值 ee 
E(D= 2 2 IPs; 
a. 


系统 中 已 坏 的 服务 台数 的 期 望 值 E(D) = N -~ E(B) - E(7)。 
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M7*/MI/N Repairable Queue System 
ZHU Yi-jun, BAO Yuan-yuan 


(Faculty of Science, Jiangsu University, Zhenjiang 212013) 

Abstract: This paper deals with the M”/M/N queue system subject to breakdowns where the server 
lifetimes are different between busy period and idle period. We derive the balance-equation of the 
steady-state. Because of the complexity of calculating the steady probability for N > 2, we program and 
derive the steady probability by using the Mathematica software. After getting the steady probability, 
some important queueing measures are obtained: the distribution of the effective servers; the moment 
generating function of the steady-state queue length; and the mean value of the queue length. 
Keywords: repairable queue; batch arrival; steady probability; moment generating function 


